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Transverse momentum dependence of Hanbury Brown-Twiss radii of pions from a 
perfectly opaque source with hydrodynamic flow 
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We investigate tlie transverse momentum dependence of pion HBT radii on the basis of a hy- 
drodynamical model. Recent experimental data show that -Rout /-Rsidc < 1, which suggests a strong 
opaqueness of the source. In addition to the opaqueness naturally caused by transverse flow, we 
introduce an extrinsic opacity by imposing restrictions on the pion emission angle. Comparing the 
HBT radii obtained from the normal Cooper-Frye prescription and the opaque emission prescription, 
we find that -Rout /-Rsido is less than unity only for small values of the transverse momentum with 
an opaque source. However, HBT radii for large values of the transverse momentum are dominated 
by the transverse flow effect and are affected less by the source opaqueness. 

O ■ PACS numbers: 25.75.Gz 

' Keywords: Relativistic Heavy Ion Collision, HBT effect, Hydrodynamical model, freeze-out 

(N : 
a : 



o 
o 



:3 



X 



I. INTRODUCTION 



I Pion intcrferomctry is one of the most promising tools for use in high energy heavy ion collision experiments designed 
J> ' to explore the states of matter under extreme conditions. As is well known, the Hanbury Brown-Twiss (HBT) effect, 
, caused by the symmetry of the wave function of identical bosons, provides information concerning the geometry of 
' the source via two-particle intensity correlation functions j^]- Recently, the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) has started to operate at extremely high energies. This opens a new frontier 
of heavy ion collision experiments. Experimental data obtained at the partial collision energy ^/s = 130A GeV and 
I preliminary data at -^s = 200 A GeV have been reported Q. One of the most interesting but puzzling results is 
' that derived from the pion HBT data Usually, the outward HBT radius -Rout is believed to be larger than the 
_ sideward HBT radius -Rsidc by an amount equal to the duration time, i.e. -R^^j ~ R^^^^ -\- Pj^{At)'^ . A hydrodynamical 
, model with a phase transition naively predicts large outward HBT radii around the RHIC energy because of the 
' I ' prolonged lifetime of the fluid due to the existence of the phase transition |^ . The obtained HBT radii in the RHIC 
O , experiments, however, reveal the surprising feature that -Rout is almost the same as (or even smaller than) -Rsidc- The 
ratio -Rout/^sido, which is proposed for a good indicator of the long emission duration [^, exhibits a slight decrease 
around unity with the pair transverse momentum. This strange result is sometimes called the "HBT puzzle" . 

Let us briefly survey the present situation regarding HBT radii. The collision process is governed by strong and 
multi-body interactions, including multiparticle production. At this time, we are far from the dynamical description 
of the entire process in terms of the fundamental theory. Though pions possess information only at their freeze-out, 
' because of the strong interaction, the two-pion correlation function provides the space-time distribution of the freeze- 
out point and the history of the space-time evolution, which is subject to the equation of state. Therefore, a dynamical 
model is indispensable for understanding the HBT data, and a hydrodynamical approach is quite suitable for this 
purpose. In addition, recent experimental data concerning an anisotropic flow (^2) obtained in the mid-rapidity 
region strongly support the validity of the hydrodynamical picture at the energies typically studied with the RHIC 
0. However, conventional hydrodynamical model analyses that reproduce single-particle spectra (and elliptic flow in 
Ref. 0) result in unsatisfactory HBT radii Space-time evolution with a smooth crossover transition equation 

of state provides a small improvement, but the resulting predictions for the HBT radii are still far from agreeing 
with the data 0- Though smaller HBT radii can be obtainedjespecially in the l ong itudinal direction, by introducing 
the chemical freeze-out in addition to the thermal freeze-out pJj] (see also Ref. |il|), this decrease is still insufficient 
to account for the experimental results. Some transport calculations indicate that the creation of dense partonic 
matter has a tendency to improve the HBT radii 13]. The problem will be the description of the hadronic stage and 
the subsequent freeze-out. A possible solution of this problem may be to introduce a modification of the freeze-out 
hypersurface in hydrodynamical models, on which the number of the emitted particle is usually carried out by using 
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the Cooper-Frye prescription |l4| with a sharp three-dimensional hypersurface. A more sophisticated treatment of 
the freeze-out combined with the hadronic transport calculation is still not sufficient to obtain predictions consistent 
with experiment . 

Considering the meaning of the HBT radii given by second order moments of the emission function T^, the 
experimentally obtained results that i?out ^ ^sidc indicate an "opaque source" |l7l IT^ . from which particles are 
emitted only on a thin surface. It has been shown that transverse flow in the hydrodynamic evolution also naturally 
causes opaque features of the HBT radii T^, and such effects are already automatically taken into account in the 
calculation. However, this effect is not sufficient to realize consistency with the experimental data. Therefore, we 
must introduce additional mechanisms that increase the source opacity.^ In the present paper, we investigate the 
HBT radii on the bases of a hydrodynamical model 0. By putting a restriction on the emission angle of pions, we 
introduce complete opaqueness of the source given by the hydrodynamical model. Comparing the HBT radii obtained 
from normal emission with those from this opaque source model, we study the transverse momentum dependence of 
the HBT radii and clarify its origin. 

In the next section, we briefly explain our model. In 3IIII a description of the opaque source model is given. Section 
IIVI is devoted to results and discussion. 



II. MODEL 



We describe the space-time evolution of hot matter created in y/s = 200 A GeV Au+ Au central collisions at the RHIC 
using the hydrodynamical model presented in Ref. We can fit single-particle distributions of charged hadrons, such 
as the pseudorapidity distribution from the PHOBOS collaboration 21], the net-proton rapidity distribution from the 
BRAHMS collaboration and the identified transverse momentum distributions from the PHENIX collaboration 
|23], by adjusting parameters in the initial matter distributions. Here, we adopt the following functional form of the 
initial net baryon number distribution: 



nB{To,r],r) = ubo <. exp 
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This form is employed in order to reproduce the flat net-proton distribution for |y| < 1 obtained by the BRAHMS 
collaboration [23|, The parameter set is listed in Table. The particle distributions are calculated using the Cooper- 
Frye prescription jl3|.^ We also take into account the resonance decay contribution, as in Ref. 0. Note that our 
calculation does not reproduce absolute numbers of kaons and (anti-)protons, because we assume the single freeze-out 
condition. Incorporating chemical freeze-out (KX. ilL . i24ij , or an additional assumption for the initial stage that leads 
to stronger transverse flow [2^, makes it possible to reproduce both the px slope and the yield of each particle species. 
But such modifications do not affect our main argument below. 



III. TWO-PION CORRELATION FUNCTION FROM AN OPAQUE SOURCE 



The two-pion intensity correlation function for a chaotic source is given by |2q | 

r(nK) 1 + l^(g'-^)l' (O^ 

1(0, fci)/(0, k2) 

where ki and k2 are the on-shell four-momenta of the two emitted pions, and q and K are the relative and average 
four-momentum defined as q — ki — k2 and K = (ki + fc2)/2, respectively. We choose I{q, K) as 

I{q,K) = J^K-da{x)e''^-y{x,K), (3) 



^ In Ref. |2Cll . HBT radii are investigated with an opaque source model inferred from quark and gluon evaporation. Their analysis given 

there, however, is based on (l+l)-dimensional longitudinal expansion and ignores transverse flow. 
^ In the present calculation, we include not only the time-like component of dtTp, as done in Ref. but also the space-like component 

of surface elements. 
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where f{x,K) is the Bose-Einstein distribution function [exTp{K ■ U{x)/T{x)) — so that l{0,k) reduces to the 

Cooper-Frye formula of a single-particle distribution |23| • Note that K is the off-shell four- momentum, and S denotes 
that the integration is carried out over the 3-dimensional freeze-out hypersurface determined by T = Tf. 

As mentioned in SjJ the opaqueness of the source represents a possible solution of the HBT puzzle. Here, we 
introduce a simple opaque source model by imposing two conditions on the emission from the freeze-out hypersurface. 

The first condition is the space-time geometrical constraint that a particle can not be emitted into the fluid; i.e., 
"inward" emission is forbidden. The assumption of local thermal equilibrium implies isotropic particle emission in 
the local rest frame of a fluid element. Therefore, inward emission naturally occurs even after a Lorentz boost to the 
center-of-mass system, though its contribution to the multiplicity is known to be small. In the Cooper-Frye freeze-out 
prescription, the number of particles absorbed from time- like surfaces can be considered as negative H^, |2^ Hfl |31| . 
We omit such emissions by introducing the step function 0{p ■ da) for p = ki and ^2 in the surface integrations. (See 
also Fig. HI) 

The second condition is a purely three-dimensional geometrical constraint that is characterized by the factor 
6{cos{(p — where (p = tan''-^ (ry/rx) is the azimuthal angle of the emission point, V' = ta.n^^{ky/kx) is the 

azimuthal emission angle of the particle (Fig. I^J, and is the step function. This constraint prohibits the emissions 
in the direction opposite to the radial flow velocity, which occur naturally by virtue of the assumption of the local 
thermal equilibrium. For example, for a particle with k^ = fcx and ky = 0, we allow emissions from space-time points 
at r,j; > 0, which corresponds to the limited range of azimuthal angles of emission points — 7r/2 < (f> < tt/2 . Hence, 
we regard the source constructed with these two conditions as a perfectly opaque source, from which particles cannot 
be emitted from the side opposite to the detector. 

Note that the above constraints should be imposed on each particle, while the two-particle correlation function 
can be expressed as the Fourier transform of the (pseudo-) single-particle distribution function. (See also Fig. [2Ib).) 
Therefore, we input the constraints on each emitted particle into both numerator and denominator of Eq. ^ as 



where Q{x,ki) = 9{ki ■ da)9[cos{4> — tpi)] and 4'i = is,in.~^{kiy/kix) for i — 1,2. Under these constraints, energy 
conservation between the fluid and emitted particles through the freeze-out process is violated. We have estimated 
the loss of energy in a system with the conditions described above to be 11%. This means that 89% of emitted energy 
quanta satisfy the conditions. The resultant multiplicity dN / drjp at rjp = 0, where rjp denotes the pseudorapidity 
rjp — 0.51n[(|p| +Pz)/{\p\ ~ Pz)], decreases from 618 to 593. This loss can be considered neghgible for the discussion 
below. Some attempts have been made to avoid this kind of violation by improving the freeze-out prescription 
but it would be a formidable task to incorporate both the constraints and the conservation low in a self-consistent 
manner, and this is beyond the scope of this paper. 



In this section, we compare the "normal" HBT radii calculated using Eqs. lO and © and the "opaque" case 
calculated using Eq. Q). We focus on i?out and i?sidc, because the purpose of this paper is to clarify the effects of 
the transverse dynamics, i.e., the flow effect and opacity effect. -Rv,ng does not reflect such effects, because its 
dependence originates mainly from rapid longitudinal expansion |32| . In the numerical evaluation of the correlation 
function, the experimental window effect was simply ignored for the better understanding of the opacity effect. The 
following approximate expressions for the HBT radii in terms of second-order moments of the source function are 



This condition can also be expressed as (?(fcT • ■''t)i where rx is the radial vector in the transverse plane. 
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IV. RESULTS AND DISCUSSIONS 
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convenient [33l |: 



iiout(fc) 




(6) 



(7) 



Here 



{A{x)){k) 



k ■ daf{x, k)A{x) 



(8) 



/ k ■ d(Tf{x, k) 



where x = x — (x) and f3± — k^/Ek. These are good approximation for both of the emission prescriptions. 

Figure O displays the results for HBT radii together with recent preliminary experimental results from the STAR 
IP and the PHENIX "35^ obtained from = 200 A GeV Au+Au colHsions. 

The results of the normal freeze-out prescription exhibit deviations similar to those seen in the 130A GeV results 
0. In particular, i?sido is smaller and decreases less steeply with increasing Kt, i?out is larger, and i?out/-Rside is much 
larger, which corresponds to a long emission duration due to the phase transition. Despite the strong restriction on the 
emission direction by Eq. Q , only a slight improvement of i?sidc and i?out is seen. In the opaque source model, -Rgido 
becomes larger at small Kt, and its A't dependence becomes slightly steeper than in the normal emission case. i?out 
becomes slightly smaller at small Kt, though no improvement can be seen beyond Kt ~ 0.5 GeV/c. Consequently, 
the ratio of i?out to i?sidc improves to a value less than unity in the smallest Kt bin. Nevertheless, in the larger 
Kt region, the opaque emission does not yield results consistent with the experimental data. These results suggest 
that the emission angle restriction does not affect the correlation function at larger Kt- After all, the flow effect 
dominates in this region. This fact can be intuitively understood by calculating r^-t correlations and (r~r^) (Fig. 0} 
and by plotting the source functions in the r^-t (Fig. [SI and r^-ry (Fig.lBJ planes. The source functions are calculated 
as 



Each function is calculated at mid-rapidity and is normalized such that its maximum value is unity. For Sxti^x, ^; ^t), 
integration over the space-time rapidity r/ = 0.5 ln[(t -f- z)/{t — z)] is carried out with \ri\ < 1 in order to obtain a clear 
emission probability distribution. For the opaque emission model, Q{x, k) is inserted into the above definitions as in 
Eq. Q. 



Figure 01 reveals a large decrease of Ar^; = y (r^^) at small fcx in the opaque emission model, as expected. This 
behavior is clearly caused by the emission angle restriction factor 9{cos{<j)—^i)), which cuts off emission from x < for 
kT — (/ctjO), as seen from the left column (fcx = 50 MeV/c) of Fig. El However, such a distinct difference disappears 
at higher fcx, because the flow causes a strong suppression and enhancement through the Boltzmann factor (right 
column). Thus, as fcx increases, the relative effect of opacity becomes weaker in the presence of transverse flow. 
Another factor, —2(3±{fxt) , in i?out exhibits strong dependence on /ct, as seen in Fig. 01 though it does not affect 
-Rout greatly. In the /ct region satisfying (Ara;)opaquo ^ (Ara;)normai, the prefactor is so small that the second 
term of Eq. Q gives only a small contribution to i?out- The decrease of {fxi) in the opaque emission model is due 
to the larger {tx), which is a natural consequence of emission only from the region Tx > 0. The leading contribution 
in this region, however, is Ai, which is 6 fm/c for kT = 200 MeV/c (Fig.[7|). Therefore, the reduction of {fxi) does 
not reduce i?out so much. From Fig. El some increase of the emissivity at the edge can be seen due to the cutoff of 
time-like surface emission. (Also, as seen from Fig.^ surface elements at the edge naturally have a time-like part.) 
This fact results in the slight increase of i?sidc f^or small fcx, because i?side can be considered the width of the source 
along the y direction. 

In the present paper, we have demonstrated that a naive opaque emission model in which we forbid emissions 
through dense media does not account for the "HBT puzzle" if collective transverse flow exists. The opaqueness 
caused by the dense matter preventing the pions from passing through such media affects the HBT radii only at small 
transverse momentum. As a result, smaller values of i?out than Rside are obtained only for small kT- Incorporating 
the opacity effect transforms the shape of the source function. Nevertheless, transverse flow dominates the source 
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function for large fcx; modification of the source function by the opacity is so slight that the fcx dependence of the 
HBT radii is still dominated by the transverse flow effect. 

As mentioned above, the present study constitutes a trial studying the modification of the source shape. One can 
consider other possibilities. For example, a viscosity correction 36] and p meson broadening "sT] have been examined. 
It also should be noted that a thermal model analysis yields acceptable agreement (SSJ. This agreement results 
from the fact that their model used in that study has a positive rx-t correlation due to the choice of the freeze-out 
hypersurface. However, that surface was simply added by hand, and is not the result of a dynamical calculation. 
Further investigation is required to solve the puzzle."* 
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FIG. 1: A schematic view of the freeze-out hypersurface (thick curve) in the r-t plane (r; = 0) and the directions (arrows) of 
emitted particles, (a) and (b) are allowed directions, but (c) is forbidden. 




k2 

FIG. 2: Emission angle constraint in x-y plane. In (a), the filled region represents the allowed region that satisfies the constraint. 
Figure (b) describes the two-particle case, where the overlapping area indicated by the darkest shading is that in which the 
constraint is satisfied for both particles. 
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FIG. 3: HBT radii for 200 A GeV Au+Au collisions as functions of the averaged transverse momentum Kt- From top to 
bottom, the open squares represent the results from our model with normal freeze-out, given by Eqs. (|5J and l^. The results 
from opaque freeze-out [Eq. (|1J] are represented by the open circles. The solid circles and squares represent the STAR and 
PHENIX experimental results, respectively. 




0.2 0.4 0.6 0.8 1 



kj [GeV/c] 



FIG. 4: rx-t correlations and source width in the outward (x) direction. 
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FIG. 5; Source function in the rx-t plane. The upper graphs are for the case of the normal freeze-out procedure, and the lower 
ones are for the opaque emission case. For each freeze-out, the small transverse momentum case (/ct = 50 MeV/c, left row) 
and large transverse momentum case (A;t = 500 MeV/c, right row) are presented. 
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FIG. 6: Source function in the Vx-Vy plane. As in the previous figure, the normal and opaque emission cases are considered for 
A;t = 50 MeV/c and 500 MeV/c, respectively. 
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FIG. 7: -Rout calculated from Eq. Q and the time duration, At. 



Tables 



TA BLE I: Parameter set for 20071 GeV Au+Au collisions at RH IC. 

£max UBO TjO Or O-r] r\D Op T[ 

6.9 GeV/fm^ 0.22 fm"^ 1.3 1.0 fm 1.4 2.6 0.92 130 MeV 



